Experimental investigation of dynamical invariants in bipartite entanglement 
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The non-conservation of entanglement, when two or more particles interact, sets it apart from other dynamical 
quantities like energy and momentum. It does not allow the interpretation of the subtle dynamics of entangle- 
ment as a flow of this quantity between the constituents of the system. Here we show that adding a third party 
to a two-particle system may lead to a conservation law that relates the quantities characterizing the bipartite 
entanglement between each of the parties and the other two. We provide an experimental demonstration of this 
idea using entangled photons, and generalize it to N-partite GHZ states. 



Introduction. Proper understanding of the production, 
quantification and evolution of quantum entanglement has 
been a major challenge of quantum physics, with direct im- 
plications on the relevance of this resource for applications in 
quantum information. A particularly important problem con- 
cerns the decay of initially entangled states under the influence 
of independent reservoirs acting on each part of the system. 
While each of these parts undergoes a typical decoherence 
process, affecting the populations and the coherences of the 
state, the dynamics of entanglement may differ considerably 
from local dynamics ill- 151. 

Usually, entanglement is not a conserved quantity. For in- 
stance, when an initially excited atom decays, releasing a pho- 
ton into a zero-temperature environment, the initial and final 
states of the atom-environment system are not entangled, but 
the atom does get entangled to the environment at intermedi- 
ate times. We show nevertheless that adding a third party to 
this two-party system leads to a conservation law involving 
quantities that measure the bipartite entanglement between 
each part of the system and the other two parts. The introduc- 
tion of this external party unveils therefore a hidden conserva- 
tion law for entanglement. We demonstrate it experimentally 
using twin photons, which have been useful tools for explor- 
ing subtle properties of quantum physics and quantum infor- 
mation 1 16]. We follow here the strategy of iflll [Till , which 
allows one to study the detailed dynamics of entangled states 
with optical interferometers. We also point out a generaliza- 
tion of this conservation law to N-partite GHZ states. 

Invariants in two-qubit dynamics. Let 5 be a qubit and R 
the reservoir with which it interacts after t = Q. The initial 
product state of the S-R system is assumed to be 



Psi?(0) =/9s(0)®Pi?(0), 



(1) 



with 
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The matrix ps is written in the basis {\g) , |e)} of the ground 
and excited states of S, and \(po) stands for the ground state 
of the reservoir R. At t = both S and R start to interact in 
such a way that the following transformation holds: 
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where p — p{t) E [0, 1] is a time-dependent parameter such 
that p{0) — 0, p{oo) = 1, and |0i) denotes the first excited 
state (orthogonal to |(/)o)) of R. The map ^ corresponds pre- 
cisely to the amplitude damping channel. For different pa- 
rameterizations p{t), the transformation ^ represents several 
physical processes such as the spontaneous emission of a pho- 
ton by a two-level atom in a zero-temperature electromagnetic 
environment, or the interaction of a two-level atom with a sin- 
gle mode of the electromagnetic field inside a cavity. Accord- 
ing to the map (O, the matrices (|2]l evolve into: 



Psip) 

pr{p) 
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(4) 



where in the expression for pr{p) represents empty rows 
and columns corresponding to the infinite remaining null ma- 
trix elements. Inspection of the reduced density matrices 
shows that the information initially contained in the system 
S is transferred to the system R. The transfer is complete at 
p = 1, when the states of S and R become exchanged. 

An important observation regarding the dynamics imposed 
by the transformation (|3) is that the mean number (TV) of total 
excitations. 



(N) = {hs{p)+nR{p)) , 



(5) 
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is conserved through the entire evolution. Thus, this restricts 
the way the populations (in the referred basis) are transferred. 
Here hg and fiji are the excitation-number operators of the 
system S and R respectively. For the specific case in which 
S represents a two level atom interacting with one mode of 
the electromagnetic field in a cavity, these operators are given 
by fis = ^(I — az) and = a)a. The conservation of 
{N) follows immediately from the expressions for {nji{p)) 
and {ns{p)), 

(nsip)) = Tr[ps{p)ns]^ Pee{l-p), (6) 

(flRip)) = Tr[PB.{p)nR] = PeeP, 
so that (N) = Pee- 

In the following, we investigate how the conservation (|5]l 
manifests itself in the evolution of the purity tt^ =Trp^ of the 
subsystem i — S,R. With the aid of Eqs. (|4]i and (|6]l it is 
straightforward to show that Tri{p) can be written as 



TT,{p)=2{h,{p)f -2{h,{p))A + l, 



(7) 



where A = 1 — \pge\ / Pee- Inverting this expression, we ob- 
tain 



(n, (p)) ^ A ± 1 ^A2_2(l-7r,(p)). (8) 



Hence we can rewrite the conservation equation (|5} as 

Pee = A±Ws(A,p)±Wfl,(A,p), 



(9) 



where we have defined Wi{A,p) = \\/ - 2(1 - i^iij))). 
As seen in Eq. (O, the appropriate choice of the sign in front 
of >Vi depends on whether {hi{p)) is greater or smaller than 
A/2. If the initial value p^e is smaller than or equal to A/2, 
then the restriction {hi{p)) < p^e implies that {ni{p)) < A/2, 
and both minus signs should be taken in (|9]l. On the other 
hand, if p^e is larger than A/2, then the signs in (|9|l depend on 
the value of p. From Eqs. (|6]l it follows that {hs{p)) ^ A/2 
whenever p^ 1 — A/2pee, whereas {n]i{p)) ^ A/2 for every 
p ^ A/2pee- In both cases the upper/lower inequality sign 
determine the ± sign that should be used. 

With the previous results we see that as a consequence of 
the conservation of (N), the purities tts{p) and ttu{p) evolve 
in such a way that the right hand side of Eq. (|9]l - or more 
generally any function of this argument - remains constant 
during the evolution. 

Entanglement conservation. We recognized that a dynam- 
ical invariance property of the map (l3) can be used to reveal 
that more complex quantities (in this case the purity of each 
subsystem) can be combined into new invariant forms. The 
purity of a system provides a different physical content from 
the excitation number, and this is the starting point for the fol- 
lowing analysis. If the quantum state ps has purity different 
from one, then there exists a larger system in a pure state, so 
that the reduced density operator corresponding to S is pre- 
cisely the mixed state ps- We will refer to the system needed 
to purify the enlarged system as M. The recognition of this 



latter system is crucial, not only for the understanding of the 
physical origin of states such as ps(0), but also because it 
throws some light into the global dynamical properties of (bi- 
partite) entanglement for systems undergoing the interaction 
modeled by (|3]l. 

Let us suppose that ps (0) results from partial tracing over 
the system M on the pure general state 

\^{Q)) = a \Ah)\e) + /? \Ah)\g) + 7 \M{)\g) + 5 |A/o)|e), 

(10) 

where |Afo) , \Mi) are two orthogonal states of M. In this 
case the elements of the initial density matrix ps (0) are p, 
\a\^ + ,Pge = I3S* + a*-f, and pgg = 
t ~ we allow system S to interact with the environment, 
according to the transformation (l3). Then, the initial tripartite 
state |*(0)) = IV-'(O)) |0o) evolves to 



' + l7l'. At 



\^{p))^ a|Afi)(yT^|e)|0o> + 
+ ,5|Afo>(yT^|e>|0o> + 
+ (/3|Afo)+7|Mi))lff)l0o 
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(11) 



Using (fTTT l. the density matrix pmsr{p) of the complete sys- 
tem may be constructed, and the reduced density matrices 
Pi [p) with i ~ M, R, S can be computed. Since M does 
not interact at all pM is constant and given by 



Pm 



(12) 



We return to Eq. (|9]l and observe that, given it was obtained 
from a property of the interaction between R and S only, it 
remains valid once the system M has been considered. More- 
over, as seen from Eq. (fTTt . the coefficients (a and S) that de- 
termine Pee are precisely the coefficients responsible for the 
entanglement between M and the rest of the system. This is 
an important observation since it relates pee directly to such 
entanglement, and hence allows us to relate pee with the (con- 
stant) purity of the system M as follows 



A ± WA2 - 2(1 - ttm) 



(13) 



A±Wm(A). 



Once more, the ± sign depends on whether pee is larger (plus 
sign) or smaller (minus sign) than A/2. Introducing this last 
expression into Eq. (|9|l leads to 



±Wa/(A) 



±Ws{A,p)±Wr{A,p). 



(14) 



As has been stated above, the appropriate choice of signs 
for each W,; is determined by the magnitude of pee relative to 
A/2,aswellasby the value of p. An inspection of all the valid 
combinations leads finally to the following cases (we omit the 
dependence on A) : 



+ Wm = Wsip) + Wb.{p), P e [0, 1], (15) 



whenever pee £ ^ , and 

- Wm ={ m{p) + Wr(j>), P e [1 - 5^] (16) 



,Ws(p)-Wfl(p), P> 2^1' 



vely, 



whenever pee > f - ■ 

Since the tripartite state remains pure during the evolution, 
each vTi may be regarded as a quantitative measure of bipartite 
entanglement between system i and the rest. Alternative! 
we may use a measure based on the Schmidt weight K,, 
which are related to vr^ according to Tii^p) ~ K^^{p). As 
a result, Eqs. dTsll-dTSIl stand as conservation relations for bi- 
partite entanglement. We conclude that, once the original S-R 
system is enlarged to include system AI needed for purifica- 
tion, the conservation law (|5]l acquires a new significance in 
terms of conservation of bipartite entanglement. The initial 
entanglement between S and M turns into entanglement be- 
tween M and R, and the entanglement transfer is described 
by Eqs. O-CSll. 

In the particular case in which the initial density matrix 
ps(0) is diagonal so that A = 1 (for example if 5 = 7 = 
in Eq. (fToli). then the conservation relations (fT5]l-(fT6]l involve 

the quantities W.,{p) = W^{l,p) ^ \J 2KdF) 
explicitly written as 



i, which are 
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(17) 



as follows from the definition of Wi and Eqs. (|6]l-(|7]i. 

The conservation relations (fT5t-(fT6t for A = 1 hold even 
when we consider system M to be in an initial state entangled 
with N qubits Sj (j — 1, 2, ...N). This allows us to find new 
global quantities that are conserved during the evolution. For 
example, the expression equivalent to Eq. (flST l for the + 1 
GHZ-type state 



is: 
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(18) 
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Experimental realization. We verified the validity of 
Eqs. JTSl l and ( fT6] l using polarization-entangled photons 
generated from spontaneous parametric down conversion 
(SPDC). The polarization entanglement is prepared with a 
two-crystal source lIlSll . The experimental setup is shown in 
figure [T] A 325nm cw He-Cd laser is used to pump two 1mm 
long type-I BBO crystals. The down-converted photons are 
spectrally filtered with lOnm bandwidth interference filters, 
and spatially filtered through 2mm detection apertures, before 
detection with single-photon counting modules. Identifying 




V,W w„ PBS 




FIG. 1 : Experimental setup. 

the polarization state of photon 2 as the system S and the po- 
larization state of photon 1 as the third party M, the source is 
set up to produce the initial state 



1$) 



a\V)j,,\V)s + P\H)j,,\H),.. 



Let us further identify the longitudinal spatial mode of pho- 
ton 2 as the reservoir R, and call the initial spatial mode |0). 
A displaced Sagnac interferometer with a nested wave-plate 
can be used to implement the following transformation on the 
polarization and spatial mode of photon 1 111 , 13i[l4ll 



1^) |0> 
\V) |0) - 



■\H] 
cos 



|0), 

9\V)\0) 



sin 6* I m ID 



(20) 



where |0) and |1) refer to different spatial modes and 6 is 
twice the angle of the half-wave plate. It has been demon- 
strated with quantum process tomography that this interfer- 
ometer implements the amplitude damping channel with fi- 
delities as high as ^ 0.95 il4ll . Identifying the polarization 
states {\H) , |V^)} with the system states {\g) , |e)}, the spa- 
tial modes {|0) , |1)} with the reservoir states {It/fo) , I0i)}^ 
and p = sin^ d, the transformation (l20t becomes equivalent 
to @, provided < 9 < 7r/2. The initial state l*)^^,^ |0)^ 
evolves to 



\m) = f3\H)^,\H)g + acose\V),,\V)s)\0)j, 



^asm9\V);,,\H) 



(21) 



after propagation through the interferometer The final state 
is equivalent to given in equation ( fTTI ) with j = S = 

0. It is important to note that partially tracing over any two 
of the three subsystems of the state (l2Tl i leads to a diagonal 
reduced density matrix. Then, the purities ttj (or equivalently, 
the Schmidt weights Kj) and each Wj {j = S, M, R) term in 
Eq. ( fTSl l can be determined by local population measurements, 
made with the detectors and 1 shown in Fig. 

By rotating the half-wave plate (HWP) in the pump beam, 
we selected different values of a and (3. For instance, we se- 
lected |ap = Pee = 0.31, 0.5, 0.73 with corresponding puri- 
ties 0.97, 0.94 and 0.96, so that the conservation laws could be 
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FIG. 2: Experimental results for pee = 0.73, A = 1. The 
coloured lines correspond to the functions Wr{p) — Ws{p) (black), 
Ws(p) + Wii{p) (red) and Ws{p) — Wb.{p) (blue). The continuous 
line represents the invariant Isr- 



applied directly. These purities were calculated from quantum 
state tomography of the initial states (after passage through 
the interferometer with 9 = 0, see below). Imperfect purity is 
probably due to spatial walk-off in the crystal and imperfect 
alignment of the interferometer. These parameters indicate 
that the experimental state is quite close the the ideal initial 
pure state. 

Projective measurements on S and M are performed us- 
ing wave-plates and polarizing beam splitters to project onto 
polarization states, while projection onto spatial modes is per- 
formed by placing detectors in mode or 1. We performed 
measurements for several values of p = sin2 6'(0<6l<7r/2) 



characterizing the amplitude damping channel in 

Figure |2] shows the theoretical curves and the experimen- 
tal data of each of the three functions that define the invari- 
ant IsB. = ^ — y^M, according to Eq. (fTSl l for the case 
Pee = 0.73 and A = 1. The black curve corresponds to 
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FIG. 3: The invariant Isr for different values of pe 



Wr{p) — Wsip), the red one to Ws{p) + Wr{p), and the 
blue one to Ws{p) — VFfl(p)- The invariant (piecewise) func- 
tion Isr is here stressed by the continuous line and corre- 
sponds, as follows from Eq. ( fT6] l (with A = 1), to one of 
the above curves depending on whether p is in the interval 



0,1 



2p, 



1 



2pee ' 



( 2pee ' 



Figure|3]shows both the theoretical and experimental curves 
of the invariant quantity Isr for three different values of pee 
(in all cases A = 1). The red and blue curves, correspond- 
ing respectively to pee = 0.31 and pee = 0.5, represent the 
invariant sum Isr = ^ + Wm = Ws{p) + Wr{p) (see Eq. 
(flSll), whereas the black curve, corresponding to pee — 0.73, 
represents the invariant Isr — \ ~ Wm given by the expres- 
sion (fT6] l. From figures |2] and |3] we see that the experimental 
data fit the theoretical curves within the precision of the mea- 
surements, thus demonstrating experimentally the conserva- 
tion relations (fT5]l-(fT6ll. 

In conclusion, we showed that by adding a third party in 
the interaction between two two-level systems, it is possible 
to obtain quantities that are invariants in the evolution of the 
entanglement. We performed an experimental demonstration 
of the process using entangled photons and verified the invari- 
ance of these quantitites. 
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